Notebook for C-SAPR exploration


In [1]:
import pyart
from matplotlib import pyplot as plt
%matplotlib inline

In [4]:
filename = '/data/collection/sgpcsaprsurI7.00/sgpcsaprsurI7.00.20150506.091302.maint.cfrad.20150506_081905.871_CSAP_v1715_SUR.nc.maint'
radar = pyart.io.read(filename)

In [5]:
max_lat = 26.5
min_lat =25.4
min_lon = -81.
max_lon = -79.5

In [8]:
print(radar.fields.keys())


[u'clutter_filtered_specific_differential_phase', u'differential_phase', u'normalized_coherent_power', u'clutter_filtered_reflectivity', u'vertical_reflectivity', u'clutter_filtered_differential_reflectivity', u'reflectivity', u'clutter_filtered_spectrum_width', u'clutter_filtered_normalized_coherent_power', u'copolar_correlation_coefficient', u'clutter_filtered_velocity', u'specific_differential_phase', u'velocity', u'clutter_filtered_differential_phase', u'clutter_filtered_copolar_correlation_coefficient', u'spectrum_width', u'differential_reflectivity']

In [9]:
#create an instance of the class using our radar
display = pyart.graph.RadarMapDisplay(radar)
#create a Matplotlib figure
f = plt.figure(figsize = [17,8])
#now we are going to do a three panel plot, resolution is a basemap parameter and determines the resolution of 
#the coastline.. here we set to intermediate or 'i' ('h' for high 'l' for low)
plt.subplot(2, 3, 1) 
display.plot_ppi_map('differential_reflectivity',
                     vmin = -7, vmax = 7,
                     resolution = 'i')
plt.subplot(2, 3, 2) 
display.plot_ppi_map('reflectivity',
                     vmin = -8, vmax = 64,
                     resolution = 'i')
plt.subplot(2, 3, 3) 
display.plot_ppi_map('velocity',
                     vmin = -15, vmax = 15,
                     resolution = 'i')

plt.subplot(2, 3, 4) 
display.plot_ppi_map('clutter_filtered_differential_reflectivity',
                     vmin = -7, vmax = 7,
                     resolution = 'i')
plt.subplot(2, 3, 5) 
display.plot_ppi_map('clutter_filtered_reflectivity',
                     vmin = -8, vmax = 64,
                     resolution = 'i')
plt.subplot(2, 3, 6) 
display.plot_ppi_map('clutter_filtered_copolar_correlation_coefficient',
                     vmin = 0.5, vmax = 1,
                     resolution = 'i')



In [ ]: